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Recent models invoking extra spacelike dimensions inhabited by (bulk) neutrinos are shown to 
have significant cosmological effects if the size of the largest extra dimension is R > 1 fm. We consider 
effects on cosmic microwave background anisotropies, big bang nucleosynthesis, deuterium and 6 Li 
photoproduction, diffuse photon backgrounds, and structure formation. The resulting constraints 
can be stronger than either bulk graviton overproduction constraints or laboratory constraints. 
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In this Letter we describe several cosmological con- 
straints on models for neutrino mass which rely on bulk 
fermions propagating in compact extra spacelike dimen- 
sions. Extra spacetime dimensions, long provided un- 
der the aegis of Kaluza-Klein (KK) and superstring the- 
ories, have played an essential role in recent attempts 
to solve fundamental problems in particle physics 
In particular, in some theories invoking n compact ex- 
tra spacelike dimensions, all Standard Model (SM) fields 
are localized on a three-dimensional surface (3-brane), 
but gravity experiences the full spacetime (bulk) [2|. 
This yields the relation Mpj = Mp +2 V n between the 
new fundamental (4 + n)-dimensional reduced Planck 
scale Mp and the four-dimensional reduced Planck scale 
Mpi = (47tGjv) -1 ^ 2 , where V n is the volume of the ad- 
ditional dimensions. If the volume of the internal space 
is sufficiently large, Mp can be much smaller than Mp\, 
giving rise to a low-scale theory of quantum gravity (e.g., 
V 2 ~ 5 x lCT 5 mm 2 gives Mp ~ lOTeV for n = 2). 

In this framework, if Mp is sufficiently small, there 
is no longer a heavy mass scale available in the theory 
to suppress neutrino masses relative to other fermion 
masses via a seesaw or similar mechanism 0]. Several 
higher-dimensional mechanisms have been developed , 
however, which can give neutrino masses and mixings 
consistent with the solar, atmospheric, and accel- 

erator neutrino experiments [8j . One widely used scheme 
posits the existence of SM-singlet fermions (neutrinos) 
which propagate in the bulk but couple via Yukawa in- 
teractions with the SM-doublet (active) neutrinos on our 
brane. This setup gives up to three light Dirac neutrino 
masses fii associated with the active neutrino flavors u e , 
Ufi, and/or v r . In addition, each bulk neutrino appears 
on our brane as a tower of massive KK modes (i.e., ster- 
ile neutrinos), and the vacuum mixing angle between an 
active neutrino and a mode with mass m mo de 3> Mi is 
r? mo d — V / 2A*i/rftmodc- The mass distribution of the 
modes depends on the geometry of the internal space. 
The simplest and most widely adopted geometry of the 
internal dimensions is that of an n-dimcnsional torus with 
radii R (1 < j < n), for which the mode masses are 
TO k = ^i/ii? + ••• + fc 2 /i? 2 , where, in bulk neutrino 
models, k = (ki, ■ ■ ■ ,k n ) is an n-tuple of whole num- 



bers, and where we assume that the bare masses of the 
bulk fermions are negligible. Several authors have found 
non-standard solutions to the neutrino anomalies in this 
framework These solutions require R^ 1 < leV 

(i?i > 0.1 /im) for the largest dimension R\, for other- 
wise theyreduce to those for a standard Dirac neutrino 
mass 8, 9J. We show how these and non-toroidally com- 
pactified models with densely distributed KK modes af- 
fect standard cosmology through their production in the 
early universe and subsequent decay. 

The incoherent production of sterile-KK neutrinos of 
mass mfc in the early universe is a nonthermal process 
governed by a Boltzmann equation 
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where /, = /; (p, t) are momentum- and time-dependent 
distribution functions, and where a is an active neutrino 
label and k, I are mode labels of a specific KK tower. 
We discuss the case in which R — Ri is the radius of 
the largest extra dimension and all other dimensions are 
small enough to have no effect on low energy neutrino 
physics. We have ignored the flavor coupling of multi- 
ple towers. The first term in Eq. is the conversion 
rate from active to sterile species and the second results 
from the decay of a mode with lifetime tj.. The latter 
arises because singlet neutrinos which mix with active 
neutrinos can decay either to SM or bulk states. On the 
brane, the partial decay width of the process Vk — > 3f 
is sin 2 6» fe G|m|/1927r 3 = G|m|juf /96tt 3 and that of the 
radiative decay Vk — > ^7 is smaller by a factor 27a/87r 
[l2| . We have also included in our calculations the contri- 
butions to Tk from visible and hadronic decays estimated 
from the partial decay widths of the Z° boson ^| . In the 
bulk, the £/-summed width of the decay Vf. — > vyhk—k' 
is ~ m\R/V2,irMp l , where hk-w is a KK graviton mode 
[l4| . The last term in Eq. (JIJ represents the decay con- 
tribution of all higher modes I > k into mode k, and 
is the appropriate collision operator. 

The conversion rate T a k — T a k(p,t) — (T/2)(P a k) to 
KK modes is the product of half the interaction rate T 
of the neutrinos with the plasma and the average prob- 
ability (P ak ) that an active neutrino v a scatters into 
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The probability depends on the matter mixing angle and 
the damping rate D = T/2 |15|: 
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Here, A k ~ m\/2p, and V — V L + V T is the full weak po- 
tential including lepton and thermal contributions. We 
assume a small lepton number of order the baryon num- 
ber (so V L <C V T ), since a larger lepton number serves 
only to enhance sterile neutrino or anti-neutrino produc- 
tion. Eq. J2J) incorporates the standard physically well- 
motivated two-neutrino active-sterile matter mixing an- 
gle and the effects of quantum damping 15]. Our con- 
straints depend on the deleterious effects of the relatively 
high modes, in which regime this formalism is identi- 
cal to that derived from direct diagonalization of the 
tower mass matrix . The finite temperature potentials 
V T from the neutrino and charged lepton backgrounds 
of the same flavor are (8v / 2G^£ , j / /3m|)(p ; y + p&) and 
(8V2 G F E u /3m%f)(pi + p{), respectively 0. 

Another important effect is the dilution of modes pop- 
ulated at temperatures T ^> 100 MeV. Disappearance of 
relativistic degrees of freedom manifests as heating of the 
plasma relative to the KK modes. We include this effect 
by following separately the complete time-temperature 
relations for the photons and modes. 

We explore the cosmological ramifications of two rep- 
resentative classes of bulk neutrino models: Class I where 
M* > 250 TeV; and Class II wherein M* ~ 1 - 10 TeV. 
Here, M™+ 2 = (27r)™Af£+ 2 . In either Class the ac- 
tive neutrinos may be coupled with one, two, or three 
bulk neutrinos. The cosmological overproduction of bulk 
graviton modes limits temperatures in this scenario to 
be less than T* ~ 10( 6 ™- 15 )/(™+ 2 ) MeV(Af»/ TeV) which 
is an upper limit on the "normalcy" temperature T* 
at which the universe must be free of bulk modes [l7j . 
Models falling into Class I correspond to large T£ (> 
^eiectroweak ~ 100 GeV) , while those in Class II have very 
low T£ (< GeV). In many Class II models effectively 
only the lightest active neutrino couples to one KK tower. 
In either Class if the heaviest active neutrino couples to 
a tower then p 3 > \/ Sm'g K w 0.057 eV [l^. We calcu- 
late here the constraints that arise in cosmologies that 
satisfy normalcy temperature requirements of their re- 
spective class. However, if the radiation dominated era 
was never significantly above the decoupling tempera- 
ture Tjydcc ~ f MeV of the active neutrinos in a very low 
reheating scenario for inflation, then early universe con- 
straints cannot be placed on either bulk graviton or bulk 
neutrino production. 

For each of these classes we have solved numerically 
Eq. (JJJ with Ck.i = for the population of the N lowest 
modes, with fully self-consistent temperature evolution 
of all relevant species. We have performed this calcula- 
tion for a single KK tower; additional towers can serve 
only to enhance cosmological limits. The height N of the 
tower is the highest mode populated at the appropriate 



T». Since our calculations begin at the highest temper- 
ature T£ permitted by graviton overproduction limits, 
any adverse cosmological effects we find will imply that 
the normalcy temperature T* must be significantly lower 
in bulk neutrino models than implied from graviton pro- 
duction alone. We have conservatively incorporated the 
effects of decays in the bulk by assuming the decay prod- 
ucts' mass-energy negligibly affects the dynamics of the 
universe. For a given momentum the number of modes 
produced per active neutrino per log-interval of tempera- 
ture T ak I H depends implicitly on the mode mass m k via 
Eq. 13, where H is the instantaneous Hubble expansion 
rate |l0|. However, as shown analytically in Ref. [10(, for 
a mode non-relativistic at the decoupling temperature 
T v dcc of the active neutrinos, the energy density (given 
by m k T ak /H integrated over InT and the active neutrino 
distribution) is independent of the mode number, despite 
the dependence of the mixing angle on the mode mass. 

This result assumes dilution is negligible and depends 
on the modes not having decayed appreciably jlOt lllj . 
Under the latter assumption, and with some simplifica- 
tions, we can extend Eq. (9) in Ref. ^3] to obtain an 
analytic estimate N Vh (BBN) - 10~ 3 {m/leV) 2 (g { Jgl k ) 
for the energy density at T v & ec in a single mode k relative 
to that in an active neutrino species. The ratio (gl/g^ k ) 
approximates dilution effects. The statistical weight in 
relativistic particles in the plasma at T^dcc is gl, and 
is glk a * t ne epoch of maximal production of mode k. 
(Roughly, this maximal production epoch is related to 
mode mass as T max ~ 133 MeV(m fc /l keV) 1 / 3 10]-) Our 
numerical calculations follow in detail the simultaneous 
production, dilution, and decay of all relevant modes of 
various energies, giving a iVt, fc (BBN) dependence on k 
which is flat modulo the effects of dilution and decay. 

Population of KK modes in the early universe leads to 
a number of unacceptable effects that provide for com- 
pelling constraints. Our calculated cosmological con- 
straints differ from those in Ref.R|, but they complement 
the supernova limits of Refs. 0, 0], and the laboratory 
constraints of Ref. 0] ■ 

Class I model constraints are given in Fig. The total 
effective number A^BBN) of neutrino flavors at the BBN 
epoch must be less than that of 4 active neutrino species, 
since otherwise the predicted and observed abundances of 
the light elements are discordant [2l| . We require the KK 
tower contribution ^ N Vk (BBN) to be less than that of 
a single active neutrino flavor. Photoproduction of deu- 
terium (D) and 6 Li due to decay of modes after big bang 
nucleosynthesis (BBN) [2^ gives another constraint. En- 
ergetic cascades dissociate 4 He into excessive amounts of 
D [23j . The increase in energy density in relativistic par- 
ticles due to mode decay prior to cosmic microwave back- 
ground (CMB) decoupling can lead to suppression of the 
second CMB acoustic peak. The current limit is that 
the effective number of neutrino flavors at decoupling is 
AT, (CMB) < 13 at 95% certainty [Hj]. Measurements to 
higher multipole moments by the Microwave Anisotropy 
Probe (MAP) (reaching A^CMB) ~ 3.9) and Planck 
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FIG. 1: Cosmological constraints on Class I bulk neutrinos. 
Photoproduction and DEBRA constrain regions between the 
dot-dashed and short-dashed lines, respectively. The CMB 
constrains parameters below the labeled BOOMERanG, 
MAP, and Planck lines and above the long-dashed line. BBN 
constrains the region between the light solid lines. Parame- 
ters must lie above the heavy solid line to be consistent with 
the inferred age of the universe. The vertical lines arise from 
the neutrino oscillation and 3 H endpoint limits on neutrino 
masses H U [H d]. Also shown is the 218 fxxa Eot-Wash 
limit on the size of two congruent large extra dimensions [20| . 



(reaching A^(CMB) ~ 3.05) surveys will be able to 
further limit the relativistic energy present at decou- 
pling , or perhaps flag the fossil relativistic energy of 
bulk modes at R ~ 0.1 fm. The increase in energy density 
due to mode decays was found by summing the energy 
injected between the neutrino and photon decoupling 
epochs. Another significant constraint comes from the 
current limits on diffuse extra-galactic background radi- 
ation (DEBRA) due to radiative decays of sterile neutri- 
nos occurring between CMB decoupling and today. The 
photon background so produced must have a total flux 
per unit solid angle dF/dQ < (lMeV/E) cm _2 sr _1 s _1 
|26l 12^ . |28[ . The expansion age of the universe, tjj > 
9 x 10 9 yr, also provides a constraint on the energy den- 
sity in KK modes. We have found that the constraints 
from distortion of the CMB spectrum and the signals 
in the solar neutrino experiments from mode decays are 
weaker than the constraints above 0, |2|| [2^ . Note that 
the arguments above do not depend on the detailed mode 
structure but rather on the existence of a high density of 
modes. 
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FIG. 2: Structure formation constraints on Class II bulk 
neutrino models with numbers of extra dimensions n = 6, 5, 4 
(!1 HDM > 0.1 to the right of the slanted lines). By assumption, 
the size of the largest extra dimension is R, and all other 
dimensions have sizes <C pm. For example, the solar neutrino 
solutions of Dvali & Smirnov and Caldwell, Mohapatra, & 
Yellin 5} lie near 0. 



Class II model effects are shown in Fig. [21 Though 
high-lying modes in these models are absent owing to low 
T* , there may remain enough energy density in low mass 
modes to comprise an appreciable hot dark matter com- 
ponent. Structure formation considerations |29j suggest 
that a hot component cannot contribute f2 HDM > 0.1. 
Contours of f2™ M = 0.1 are shown in Fig. |3 for these 
models with n = 6,5,4 extra dimensions. Note that 
some recent models for solar neutrino oscillations fall 
in a parameter range which could give an appreciable 
qhdm Whether this can constitute a true constraint 
depends on the precise relation between T* and M* in 
these models [3(|. At present all Class I and II mod- 
els 0, S 0] can escape elimination by invoking a suf- 
ficiently low (< 20MeV) re-heating temperature T r for 
inflation. 

The spectrum of low-lying modes in Class II mod- 
els could give a viable dark matter candidate if T r is 
low. Low T r results in suppressed production of high- 
mass modes that provide the closure, BBN, and decay 
constraints. For low enough all modes produced be- 
low T r survive until today and escape decay constraints. 
The possibility of a realistic dark matter component 
from the KK modes is finely tuned. For instance, for 
R- 1 = 4 x 10~ 8 MeV, and fr = 10~ 5 eV, ft„ 3 - 0.1 for 
T r - lGeV, but Q„ s - 0.2 for T r - 1.3 GeV. Albeit 
finely tuned, the latter case, for which f2™ M < 0.1, is an 
interesting dark matter candidate, comprising a mixture 
of hot, warm (~ keV), and cold (~ MeV) components. 

Modifications to Class I and II models may allow cir- 
cumvention of our constraints. A stronger dependence 
of the mixing angle on M,/m mo d e would ensure that the 
population of the modes would fall with increasing mass. 
An alternate dependence of mode lifetime on m mo de and 
Hi could eliminate some or all of the constraints. There 
could exist multiple additional (possibly fat) branes in 
the bulk, devoid of energy density and parallel to our 
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own, onto which modes decay preferentially |l7j. If the 
re- heating temperature T r of inflation is near T vt ^ cc , no 
KK modes will be populated in the radiation dominated 
era, and therefore the constraints presented here do not 
apply. Some population of KK modes can occur during 
reheating, or through resonant production if there is a 
large lepton number, but we do not explore these possi- 
bilities here. Also, the internal dimensions need not be 
toroidally compactified [3l|. A space which has a KK 
mode decomposition with a sufficiently low mode den- 
sity (or a gap) could evade cosmological constraints. Bad 
effects of higher /ii can be removed by restricting the 
KK towers to those built on low [ii (< 1CU 4 eV) active 
neutrinos — this is what is (or should be) done in Class 
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